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OBJECTIVES 

To  carry  out  fundamental  and  wide  ranging  research  investigations  involving  the  nonlinear 
wave  propagation  which  arise  in  physically  significant  systems  with  emphasis  on  nonlinear  optics. 
The  modeling  and  computational  studies  of  wave  phenomena  in  nonlinear  optics  include  nonlinear 
wave  propagation  in  photonic  lattices,  ultrashort  pulse  dynamics  in  mode-locked  lasers,  dynamics 
of  dark  solitons  and  analysis  of  dispersive  shock  waves. 

STATUS  OF  EFFORT 

The  Pi’s  research  program  in  nonlinear  wave  propagation  with  emphasis  in  nonlinear  optics  is 
extremely  active  and  covers  a  number  of  research  areas.  There  have  been  a  number  of  important 
research  contributions  carried  out  as  part  of  the  effort .  During  the  period  15  April,  2012-14  April, 
2015,  fourteen  papers  were  published  in  refereed  journals.  In  addition,  two  refereed  conference 
proceeding  were  published,  and  twenty  one  invited  lectures  were  given.  The  PI  also  received  a 
number  of  honors  during  this  period. 

The  key  results  and  research  directions  are  described  below  in  the  section  on  accomplish¬ 
ments/new  findings.  Full  details  can  be  found  in  our  research  papers  which  are  listed  at  the  end  of 
this  report. 
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The  lattice  nonlinear  Schrodinger  equation,  which  is  derived  from  Maxwell’s  equations  under 
the  paraxial  approximation,  provides  the  starting  point  of  the  analysis  associated  with  lattice  pho¬ 
tonics.  The  linear  index  of  refraction  is  associated  with  the  lattice  background;  the  nonlinearity 
is  derived  from  the  inclusion  of  cubic  nonlinearity.  Often  the  background  potential  is  periodic, 
but  not  always.  In  our  research  we  have  also  studied  quasi-crystal  lattices  and  lattices  with  de¬ 
fects;  we  have  carried  detailed  computational  studies  investigating  stability  properties,  boundaries 
of  associated  bandgaps  with  solitons  and  vortices. 

When  the  Schrodiunger  potential  is  taken  large,  one  can  employ  an  approximation,  called  the 
tight  binding  limit.  Within  this  approximation  key  nonlinear  discrete  systems  are  derived.  In  a 
regime  near  ‘Dirac  points’,  associated  with  special  locations  in  the  underlying  linear  dispersion  re¬ 
lation,  which  exhibit  conical  singular  structure,  new  phenomena  are  found.  Some  of  the  properties 
can  be  described  in  terms  of  the  continuous  limit  where  special  (nonlinear  Dirac  type)  systems  are 
obtained.  When  the  background  lattice  is  weak  different  types  of  Dirac  systems  result. 

Our  first  research  efforts  describe  the  wave  dynamics  in  situations  where  boundary  effects  were 
not  important.  Subsequently  we  investigated  cases  when  boundary  edges  and  effects  are  important. 
In  the  latter  case  novel  dynamics  are  fund  to  occur.  In  such  situations  sometimes  unidirectional 
waves  associated  with  Floquet  topological  insulators  can  be  found.  Such  waves  can  propagate  over 
very  long  distances  with  no  backscatter.  This  situation  is  under  current  careful  study. 

Besides  studying  the  dynamics  and  propagation  of  edge  wave  dynamics  in  honeycomb  lattices 
we  have  also  analyzed  the  propagation  of  multidimensional  nonlinear  surface  plasmon  waves  in 
dielectric/metal  interfaces. 

Novel  continuous  and  discrete  systems  closely  related  to  the  classical  nonlinear  Schrdinger 
equation  and  which  exhibit  parity  time  symmetry  were  solved  exactly  via  the  inverse  scattering 
transform.  They  have  have  the  property  that  the  nonlinear  induced  potential  simultaneously  de¬ 
pends  on  current  and  reflected  spatial  locations. 

The  properties  of  mode-locked  lasers  which  are  used  to  create  ultrashort  pulses  were  analyzed 


2 


under  conditions  where  one  has  higher  order  perturbative  terms  such  as  third  order  and  nonlinear 
dispersion.  The  basic  mathematical  model  employs  gain  and  filtering  terms  saturated  by  energy  and 
a  loss  term  saturated  by  power;  this  is  referred  to  as  the  power  energy  saturation  (PES)  equation. 
Since  dark  solitons  in  mode  locked  lasers  have  been  observed,  we  also  used  this  model  to  study 
their  mode  locking  characteristics. 

Motivated  by  experiments  at  the  University  of  Colorado,  dispersive  shock  waves  (DSWs)  were 
also  investigated.  DSWs  are  found  in  many  applications  including  nonlinear  optics,  fluid  dynamics 
and  Bose  Einstein  condensation.  We  showed  that  for  certain  fundamental  nonlinear  wave  equations 
with  multi-step  initial  conditions,  even  though  the  data  breaks  up  in  to  numerous  DSWs  in  an 
intermediate  long  time  limit,  eventually  the  solution  tends  to  one  DSW. 
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ACCOMPLISHMENTS/NEW  FINDINGS 


Nonlinear  photonic  lattices 

Nonlinear  light  wave  propagation  in  periodic  optical  photonic  lattices  is  an  exciting  and  active 
area  of  research.  This  is  partially  due  to  the  realization  that  photonic  lattices  can  be  constructed  on 
extremely  small  scales  of  only  a  few  microns  in  size.  Hence  they  allow  the  possibility  of  careful 
manipulation  and  navigation  of  lightwaves  on  small  scales.  Localized  nonlinear  optical  pulses 
have  been  studied  when  they  occur  on  one  and  two  dimensional  backgrounds.  In  one  dimension 
the  periodic  backgrounds  are  often  referred  to  as  waveguide  arrays 

The  periodic  backgrounds  can  be  either  fabricated  mechanically  such  as  those  comprised  of 
AlGaAs  materials  (cf.  [1]  where  one  dimensional  arrays  are  considered)  or  all-optically  using 
photo-refractive  materials  as  done  in  [2,  3].  These  latter  photonic  structures  can  be  constructed  via 
interference  of  two  or  more  plane  waves. 

Waveguide  arrays 

One  dimensional  coupled  optical  waveguides  is  a  setting  that  represents  a  convenient  labora¬ 
tory  for  experimental  observations.  By  etching  the  surface  of  the  optical  material  one  can  form 
a  periodic  structure  which  is  referred  to  here  as  a  waveguide  array.  Typically,  such  a  waveguide 
array  is  composed  of  numerous  single-mode  individual  waveguides  each  being  a  few  microns  in 
width  and  a  few  millimeters  in  length.  Such  small  scale  structures  could  be  embedded  in  a  large 
scale  environment  and  potentially  can  be  used  to  guide  light  in  a  controllable  manner. 

In  general  an  input  optical  beam  will  diffract  in  the  lateral  or  horizontal  direction  unless  it  is 
balanced  by  nonlinearity.  The  first  theoretical  prediction  of  nonlinear  wave  propagation  in  an  opti¬ 
cal  waveguide  array  which  possessed  localized  discrete  waves;  i.e.  discrete  solitons,  was  reported 
in  [4].  The  wave  phenomena  satisfies  the  well  known  discrete  nonlinear  Schrodinger  (DNLS) 
equation 

dE 

i —Qz  t~  C(En+ 1  +  En_i )  +  n\En\2 En  =  0  , 

where  En  is  the  field  strength  at  site  n,  C  is  a  coupling  constant  which  depends  on  adjacent  waveg- 
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uides,  k  is  a  constant  describing  the  nonlinear  index  change  and  z  is  the  direction  of  propagation. 

In  1998  self-trapping  of  light  in  a  discrete  nonlinear  waveguide  array  was  experimentally  ob¬ 
served  by  a  group  of  researchers  at  the  Weizmann  Institute,  Israel  [1].  In  their  experiment,  an 
array  of  ridge  waveguides  were  etched  onto  an  optical  substrate.  In  the  experiments,  at  low  power, 
linear  behavior  is  observed;  see  Figure  1-top.  The  light  is  found  to  spread  among  nearly  all  the 
waveguides  hence  experiencing  discrete  diffraction.  However,  at  sufficiently  high  power,  the  beam 
self-traps  to  form  a  localized  state  (a  soliton)  in  the  center  waveguides  see  Figure  1  -bottom. 


70  w 

320  W 

500  W 


Figure  1:  Images  of  the  output  for  different  powers.  Top  figure  corresponds  to  low  peak  power; 
linear  features  are  demonstrated.  Middle  inset  is  taken  at  moderate  peak  power  where  it  is  seen 
that  the  pattern  is  shrinking.  Bottom  figure  corresponds  to  high  peak  power  At  sufficiently  high 
power,  a  localized  structure,  i.e.  a  discrete  soliton,  is  formed. 

Subsequently  many  theoretical  studies  of  discrete  solitons  in  waveguide  arrays  which  found 
switching,  steering  and  other  interesting  properties. 

This  motivated  many  interesting  properties  of  nonlinear  lattices  and  discrete  solitons  including 
switching,  steering  and  other  interesting  properties.  Discrete  systems  exhibit  phenmomena  that  are 
absent  in  continuous  media  such  as  the  possibility  of  producing  anomalous  diffraction  [5].  Hence, 
in  discrete  systems  self-focusing  and  defocusing  processes  can  be  achieved  in  the  same  medium 
(structure)  and  wavelength.  This  also  leads  to  the  possibility  of  observing  discrete  dark  solitons  in 
self- focusing  Kerr  media  [6].  The  experimental  observations  of  discrete  solitons  [1]  and  diffraction 
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management  [5]  also  motivated  interest  in  discrete  solitons  in  nonlinear  lattices. 

In  [7]  we  considered  the  possibility  of  moving  discrete  solitons;  we  found  that  as  discretiza¬ 
tion  effects  become  important  the  solitons  slowed  down.  This  work  motivated  our  studies  of  more 
complicated  situations  in  which  laser  beams  can  propagate  in  a  discrete  medium  with  alternating 
diffraction.  A  natural  question  which  arises  is:  What  happens  if  we  launch  a  high  power  laser 
beam  into  a  diffraction  managed  waveguide  array  [8]?  We  found  modes  which  “breathe”  peri¬ 
odically  upon  propagation.  During  propagation  the  peak  amplitude,  FWHM,  and  phase  of  the 
beam  evolve.  For  example,  during  each  map  period,  the  FWHM  of  the  beam  oscillates  between  a 
nonzero  minimum  and  a  maximum  value  while  maintaining  conservation  of  power. 

By  designing  the  diffraction  properties  of  a  linear  waveguide  array,  additional  new  phenomena 
are  found  to  occur.  In  [9]  we  showed  that  by  appropriately  tailoring  the  diffraction  properties  of 
a  waveguide  array,  the  interaction  between  discrete  solitons  can  be  altered  to  achieve  better  tun- 
ability  and  control  over  the  collision  outcome.  For  example,  by  colliding  vector  discrete  solitons, 
remarkably  the  interaction  picture  involves  beam  shaping,  fusion  and  fission.  Interestingly  it  was 
also  shown  that  discrete  solitons  in  two-dimensional  networks  of  nonlinear  waveguides  can  be  used 
to  realize  intelligent  functional  operations  cf.  [10]. 

Two  dimensional  periodic  nonlinear  photonic  lattices 

In  many  photonics  applications  a  nonlinear  Schrodinger  equation  with  an  external  potential,  or 
the  lattice  NLS  equation,  is  the  governing  equation.  It  can  be  derived  from  Maxwell’s  equations, 
and  in  dimensionless  form  it  is  given  by 

iuz  +  A u  —  V (r )u  +  \u\2u  =  0  .  (1) 

Here  u( r,  z)  corresponds  to  the  slowly-varying  complex  amplitude  of  the  electric  field  in  the  r  = 
(x,  y)  plane  propagating  along  the  z  direction,  A u  corresponds  to  diffraction,  V (r)  is  an  optical 
potential  representing  the  linear  index  of  refraction  and  the  cubic  nonlinear  term  originates  with 
the  nonlinear- Kerr  change  of  the  refractive  index  often  referred  to  as  the  yA  nonlinear  index. 

The  potential  is  assumed  uniform  along  the  propagation  direction  z,  although  non-uniform 
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potentials  in  z  could  also  be  considered  in  a  manner  analogous  to  diffraction  or  dispersion  man¬ 
agement,  and  more  recently  has  been  employed  in  a  novel  way  in  honeycomb  lattices  [11]. 

Most  studies  consider  V (r)  to  be  periodic  in  r  =  (x,  y ),  though  some  work  on  lattice  defects 
have  been  carried  out  cf.  [12].  When  the  lattice,  V ( x ,  y )  is  given  by 


v(x,y)  =  ^\J2ei 


N- 1 


(kxx+kyy)  1 2 


N2 


(2) 


where  kx  =  kcos(2im/N),  ky  =  k  sin(27rn/ZV),  k  are  constant,  it  is  periodic  with  iV  =  4;  this  is 
sometimes  called  a  square  lattice.  Aperiodic  lattices  or  quasi-crystals  arise  when  N  =  5,7. 

In  [13]  solitons  were  observed  in  “irregular  ”  quasi-crystal  lattices  (e.g.  N=5)  and  in  [14]  a 
detailed  computational  study  carried  out  for  this  potential  and  for  other  potentials  with  defects. 
In  [15]  we  studied  the  stability  properties  and  the  boundaries  of  the  bandgaps  associated  with 
the  above  potential  for  the  lowest  band  of  the  dispersion  surface;  and  in  [16]  we  computationally 
investigated  the  existence  and  stability  of  vortex  solitons. 

Next  we  consider  certain  two-dimensional  periodic  optical  external  potentials: 

(a) 


y  / 

C-*  -A 

_ JN+Vj 

H  So 

Vl 

's-« 

v2 

X 

-> 


Figure  2:  Left:  simple  square  lattice;  right:  non-simple  honeycomb  (HC)  lattice 


V  (r  +  mvi  +  nv2)  =  V  (r) 

with  m,n  €  Z,  separated  into  two  categories:  simple  and  non-simple  honeycomb  (HC)  lattices. 
In  a  simple  lattice  all  sites  can  be  constructed  from  one  site,  whereas  non-simple  HC  lattices  are 
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constructed  from  more  than  one  site;  Figure  2  illustrates  a  simple  square  lattice  and  a  non-simple 
honeycomb  lattice. 

In  the  case  of  a  simple  square  lattice,  the  analog  experiment  described  above  in  the  context 
of  one  dimensional  waveguides,  was  explored.  At  low  power  an  input  beam  diffracts  whereas  at 
high  power  a  2d  lattice  soliton  is  observed  [17];  see:  Figure  3-  experiment,  and  Figure  4  where 
computational  results  were  obtained  by  solving  equation  (1).  Subsequently  there  have  been  many 
studies  which  have  investigated  a  variety  of  localized  soliton  solutions  in  two-dimensional  simple 
lattices  cf.  [18,  19,  20,  21,  22] 

I 

Figure  3:  Experiment:  Left:  Low  Power;  Right:  High  Power 


Ligure  4:  Typical  2d  lattice  soliton:  Computational 

In  what  follows  next,  we  describe  some  of  our  recent  work  involving  two  dimensional  optical 
lattices. 

The  underlying  linear  problem,  is  obtained  from  equation  (1)  by  assuming  u( r,  z)  =  (^(r)e~i/'’z 
with  \ip\  <C  1,;  this  yields  the  linear  Schrodinger  equation  with  periodic  coefficients: 

(V2  -  V (r))</?  =  -w  (3) 
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where  V (r)  is  a  2-d  periodic  potential  with  lattice  vectors:  V| ,  v2  on  the  lattice  P  =  {mvi  +  nv2  : 
m,  n  G  Z}. 

Bloch  theory  implies  that  solutions  p  of  equation  (3),  called  Bloch  modes,  have  the  form  p  = 
ip( r;  k)  =  elk  rU (r;  k)  where  k  is  real,  U (r;  k)  is  periodic  in  r  with  the  same  periodicity  as  V (r) 
and  <p(-,k),ji  =  /i(k)  are  periodic  in  k  with  ‘dual’  lattice  vectors:  ki,k2  where  vm  ■  kn  = 
2n8mn',  //  =  /i(  k)  is  called  the  dispersion  relation. 

Due  to  periodicity  in  k  we  can  represent  tp( r;  k)  as  <p( r;  k)  =  £  0v(r)e*k  v  where  the  is 

V 

over  all  values  in  P;  0v(r)  is  termed  a  ‘Wannier  function’.  Using  the  physical  space  periodicity  it 
follows  that  0v(r)  =  0o(r  —  v). 

In  general  one  cannot  find  explicit  analytical  solutions  for  these  Wannier  functions.  But  when 
the  potential  is  large,  i.e.  \V\  3>  1,  which  is  termed  the  tight  binding  approximation  (TBA)  and 
is  well  studied  in  the  literature,  one  can  approximate  the  Wannier  functions.  For  simple  lattices 
Wannier  functions  are  approximated  by  “orbitals”  via: 

(V2-Uo(r))0(r)  =  -U0(r) 

where  Vo(r)  is  the  potential  in  the  fundamental  cell,  constant  outside;  (f>( r),  E  are  the  corresponding 
eigenfunction  (orbital)  and  eigenvalue.  The  approximation  involves  solving  for  0(r),  E  in  the  limit 
\V\  1;  it  is  often  useful  to  consider  Vo(r)  ~  D0(a2(x 2  +  c2y 2))  locally  near  the  minima;  in  this 

case  at  the  origin. 

In  [23,  24]  we  developed  the  necessary  mathematical  foundations  to  study  nonlinear  waves  in 
simple  and  non-simple  HC  periodic  lattices. 

Simple  lattices 

For  simple  lattices  with  one  dispersion  branch  we  can  find  discrete  evolution  equations  for  the 
envelope  associated  with  the  lattice  NLS  equation  (1)  by  looking  for  solutions  of  the  form 

u  ~  ^  av(Z)(j>( r  —  v)e*k'v 
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where  Z  =  ez,  \e\  <C  1  and  av(Z)  represents  the  slowly  varying  Bloch  wave  envelope  at  the 
site  Sv.  Substituting  this  approximation  into  the  lattice  NLS  equation  (1),  multiplying  by  <f>( r  — 
p)e_*kp,  integrating  over  all  space  after  some  manipulation  yields  a  discrete  NLS  equation  at 
general  values  of  k 


.  ClCl-p  V  A  7'k  v  ,  I  1 2  r» 

/  yQp+yCye  T  CT(^|fl-p|  Op  0 

<v> 

where  (j)  is  an  orbital  (see  above),  g  =  J[0(r)]4dr  and  <  v  >  means  we  take  nearest  neighbors 
to  point  p.  This  is  a  generalization  of  the  well-known  one  dimensional  lattice  problem  [4],  but 
extended  to  two  dimensional  square  and  triangular  lattices  [23].  These  results  can  be  extended  to 
cases  where  there  are  more  than  one  dispersion  relation  band. 

As  a  further  limit,  when  the  envelope  ov  varies  slowly  over  v  with  a  scale  R  =  z/r,  \v\  <C  1; 
then 

i/j  ~  ^  av(Z)(f>( r  —  v)e*k’v  ~  ^  a(R,  Z)<fi(r  —  v)e*k’v 

V  V 

The  continuous  limit  yields  a  continuous  NLS  equation  for  each  k  point  in  the  Brillouin  zone 

da  -  ~  v2  A  -  ~  2 

IE—  +  JA/J  ■  Va  +  —  2_^  Orn,ngom,na  +  ag\a\  a  =  0. 

m,n= 1 

where  dm  =  ,  Bm  =  ^|-  and  V/i  plays  the  role  of  the  group  velocity.  For  different  locations 

in  a  Brillouin  zone  (periodic  cell  in  k  space)  the  spatial  operator  in  the  above  NLS  equation  can 
be  elliptic,  hyperbolic  or  parabolic.  We  remark  that  further  reduction  is  possible  by  going  into  a 
moving  frame  with  the  group  velocity  and  introducing  a  slower  evolution  scale  Z  =  vZ ;  maximal 
balance  is  obtained  when  £  =  0(u)  =  O(crg). 

Non-simple  honeycomb  lattices 

Next  we  consider  background  honeycomb  (HC)  lattices  such  as  that  depicted  in  Figure  2-right; 
this  was  considered  in  photonic  lattices  by  Segev’s  group  cf.  [25,  26].  A  typical  intensity  plot  of 
a  honeycomb  lattice  is  given  in  Figure  5  where  the  local  minima  are  in  the  blue  regions,  which  in 
turn  form  an  hexagonal  lattice. 
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Figure  5:  Typical  physical  space  intensity  plot  of  a  honeycomb  lattice 


Figure  6:  First  two  bands  of  the  dispersion  relation  /i(k)  associated  with  a  typical  honeycomb 
lattice;  the  bands  touch  at  Dirac  points 

Similarly  a  typical  HC  dispersion  relation  /i(k)  for  the  lowest  two  bands  is  given  in  Figure  6.  It 
is  seen  that  the  bands  touch  (cf.  [27,  28])  at  special  points,  called  Dirac  points  which  in  turn  form 
the  edges  of  a  HC  lattice. 

In  the  carbon  based  material  graphene  Dirac  points  play  an  important  role  in  terms  of  the  novel 
phenomena  found  such  as  extreme  conductance  and  strength.  Similarly,  in  optical  applications 
new  phenomena  are  found.  In  optics  this  research  direction  is  sometimes  referred  to  as  photonic 
graphene  or  optical  graphene.  This  is  an  artificial  form  of  its  carbon  based  analogue;  it  allows  one 
to  investigate  phenomena  which  would  be  very  difficult  to  do  in  its  original  carbon  form. 

In  our  research  [29,  30]  we  have  found  novel  discrete  equations  associated  with  HC  lattices  via 
the  tight  binding  approximation.  The  solution  of  the  lattice  NLS  equation  is  assumed  to  be  of  the 
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form 


u  ~  ^2  av(Z)(j)A( r  -  v)e*k'v  +  ^  6v(Z)0s(r  -  v)elk'v 

V  V 

where  Z  =  ez,  |e|  <<  1  and  sum  v  takes  all  values  m,  n  €  Z  associated  with  A,  5  (minima  in 
each  cell)  sites  respectively  and  the  orbitals  <t>A,  On  satisfy 

(V2  -  V)(r))  4>j(y)  =  -Ej<f>j( r);  j  =  A,  B 

Substituting  u  into  the  lattice  NLS  equation  (1),  multiplying  by  <f>j(r  —  p)e_*k  p;  j  =  A,  B 
and  integrating  we  find  after  some  manipulation  and  normalization  (with  maximal  balance:  slow 
evolution,  dispersion,  nonlinearity)  the  following  Discrete  HC  system,  at  a  general  location  k: 


i  i*  +  C  bp  +  a\ap\2ap  =  0 

QjZj 

(4) 

(J  h 

+  C+av  +  a\bp\2bp  =  0 

(5) 

C~bp  =  bp  +  p&p_vie-’kvi  +  P6p„V2e-ikV2 
C+av  =  av  +  pap+vie*k'vi  +  pav+V2etk'V2 

where  a  depends  on  the  original  coefficient  of  nonlinearity  and  integrals  over  the  orbitals,  Z  = 
ez,and  p  is  deformation  parameter  which  depends  on  given  HC  lattice  parameters;  when  p  =  1  we 
have  perfect  hexagonal  HC  lattice.  The  effect  of  p  in  physical  space  is  to  squeeze  the  bottom/top 
hexagonal  comers  of  the  HC  lattice  when  p  <  1  and  push  them  apart  when  p  >  1;  when  p  <  1/2 
a  gap  opens  at  the  Dirac  points. 

The  above  system  is  a  discrete  evolution  system  which  governs  the  evolution  of  Bloch  wave 
envelopes  at  a  general  location  k.  For  different  values  of  k,  we  find  that  the  evolution  of  the 
envelopes  can  be  very  different.  We  can  consider  a  further  limit  like  we  did  for  simple  lattices;  i.e. 
dp+v  ~  (1  +  w  ■  V  H - )a;  namely  the  continuous  limit  where  the  envelopes  can  be  written  in  the 
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form  a(R,  Z)  and  6(R,  Z)  where  R  =  ur  represents  the  spatial  coordinate  of  the  slowly  varying 
envelope.  If  we  assume  the  wave  is  propagating  with  k  near  a  Dirac  point  we  find  the  following 
continuous  system  of  equations  which  is  quite  different  from  an  NLS  type  equation 

da  ~  2 

i  — — b  O-b  +  a  Jana  =  0 
oz 

db  ~ 

i-^ - d+a  +  au\b\2b  =  0  (6) 

oz 

where  d±  =  dx±i(dy,  ay  =  ^  and  (  =  \J ip2^1  ■  If  by  =  0,  find  a  linear  system  which  is  the 
familier  2  dimensional  wave  equation;  otherwise  this  is  termed  a  NL  Dirac  system. 

In  Figure  7  is  a  typical  comparison  between  the  governing  lattice  NLS  equation  and  the  above 
reduced  system  with  p  —  1.  In  Figure  7  top  panel  a-c:  the  intensity  of  the  solution  is  plotted 
while  in  the  lower  panel  d-f:  only  the  a  component  of  the  continuous  Dirac  system  (6)  is  given- 
the  b  component  looks  similar  to  a.  The  initial  condition  is  a  unit  gaussian  for  a  and  zero  for  the 
b  component.  We  see  bright  rings  emanating  from  the  center;  this  is  termed  conical  diffraction. 
Interestingly,  similar  conical  structure  occurs  in  the  linear  system-i.e.  when  a  =  0.  We  note  that 
when  the  initial  condition  for  both  a,  b  are  unit  gaussians  we  observe  conical  diffraction  with  a 
‘notch’  [29,  30], 

The  main  difference  between  the  linear  and  nonlinear  systems  is  the  phase  structure  which 
becomes  increasingly  complex  in  the  nonlinear  case.  Conical  diffraction  was  observed  experimen¬ 
tally  [25].  It  is  interesting  to  note  that  remarkable  complex  phase  structures  were  also  recently 
observed  [31]. 

In  [32]  direct  computation  on  the  lattice  nonlinear  Schrodinger  equation  (1)  showed  that  when 
the  nonlinear  coefficient  was  significantly  increased  one  obtains  triangular  diffraction  as  opposed 
to  conical  diffraction.  Using  our  discrete  system  we  showed  [33]  that  triangular  diffraction  was 
present.  Figure  8  describes  the  situation.  We  see  that  for  small  nonlinear  coefficient  the  resulting 
pattern  is  nearly  conical.  But  as  the  nonlinear  coefficient  is  increased  to  unity  a  clear  triangular 
structure  emerges.  When  the  nonlinear  coefficient  takes  the  opposite  sign  the  triangular  diffraction 
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Figure  7:  An  initial  ‘spot’  evolves  into  bright  rings;  the  top  panel  is  from  direct  integration  of  the 
lattice  NLS  equation  (1).  The  bottom  panel  is  the  a  component  of  the  above  continuous  Dirac 
system  (6)  with  p  =  1. 

changes  from  a  triangle  pointing  upward  to  one  pointing  downward. 


Figure  8:  Effective  NL  coefficient  in  the  discrete  system  with  p  =  1,  a  =  (a)  0  (b)  0.1  (c)  1  (d)  -1 

We  find  that  when  we  consider  the  continuous  limit  one  needs  to  keep  higher  order  terms  in 
order  to  see  the  triangular  structure.  The  relevant  higher  order  continuous  equations  with  p  =  1 
are  given  [33]  by 


idz  a  +  (<9_  —  ^  A+)6  +  ov|a|2a  =  0 
idz  b  +  (d+  -  ^  A_)a  +  du\b\2b  =  0 

where  d±  =  ±dx  +  idy,  A±  =  dxx  +  (hr /‘i  ±  2ic)xdy  and  the  corresponding  figures  demonstrat¬ 
ing  triangular  diffraction  are  given  below  in  Fig.  9.  If  the  higher  order  terms  A±  are  not  included 
in  the  above  equations  the  resulting  diffraction  patterns  are  conical,  not  triangular. 

Interestingly,  when  we  consider  initial  conditions  for  the  discrete  system  which  are  on  the  scale 
of  the  discretization,  then  in  the  linear  discrete  problem  with  a  =  0  we  see  spots  at  the  hexagonal 
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Figure  9:  The  nonlinear  coefficient  is  av  —  (a)  0  (b)  1  (c)  5  (d)  -5 


Figure  10:  Circular  rings  become  elliptic  when  honeycomb  lattices  are  deformed,  (a)  p  =  0.8, 
(b)  p  =  0.6. 

corners,  whereas  in  the  nonlinear  case  the  spots  are  on  triangular  corners  with  one  corner  spot 
having  larger  intensity,  upward  spot  being  of  larger  intensity  for  one  sign  of  nonlinearity  and  visa- 
versa  for  the  other  sign  of  nonlinearity  [33]. 

When  the  deformation  parameter  p  /  1  then  equation  (6)  admits  an  elliptical  diffraction  pattern 
which  can  be  seen  below  in  Figure  10.  If  we  denote  (3  =  2p  —  l,  then  for  j3  =  0(1)  and  positive 
we  have  elliptical  diffraction,  as  mentioned  above.  When  (3  =  0(1)  and  negative,  a  gap  opens  up 
in  the  spectrum  and  we  find  NLS  type  equations  which  can  contain  solitons  (gap  solitons).  When 
p  ~  1/2  and  \(3\  <§C  1  we  find  "straight  line’  diffraction.  In  fact  we  find  a  class  of  reduced  equations 
in  a  regime  where  (3  =  0(v2)  [34,  35].  One  example  of  a  member  of  this  class  of  equations  is  the 
following 


C  3 


dg(dzF -aJ\F\2F) +d^F  =  0  NLSKZ 

where  6  represents  a  right  or  left  going  characteristic  variable  obtained  from  the  leading  order 
one  dimensional  wave  equation,  obtained  when  p  =  1/2.  We  refer  to  this  equation  as  the  NL¬ 
SKZ  equation  because  of  the  analogy  with  the  so-called  classical  KZ  equation  after  Khokhlov  & 
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Zabolotskaya,  Sov.  Phys.  -Acoust.,  1969  given  by 


dg(dtu  +  udgu )  +  7  dyU  =  0  KZ 

From  a  initial  conditions  for  a,  b  of  guassian  type  we  find  that  there  are  two  small  parabolic 
structures  emanating  from  the  center.  The  asymptotic  theory  is  found  to  provide  a  good  approxi¬ 
mation  to  the  underlying  complex  situation.  More  details  can  be  found  in  our  papers  [34,  35]. 

Recently  it  has  been  shown  [11]  that  semi-infinite  HC  lattices  with  a  periodic  helical  structure 
written  along  the  propagating  direction  can  support  linear  traveling  edge  modes.  The  written  heli¬ 
cal  HC  structure  has  the  effect  of  introducing  a  pseudo-magnetic  field.  In  [1 1]  it  is  found  that  these 
linear  edge  modes  can  move  along  and  follow  the  boundary  even  when  there  are  corners.  In  pre¬ 
liminary  work  we  have  developed  an  asymptotic  theory,  assuming  a  relatively  fast  helical  change, 
to  describe  this  situation.  We  first  studied  the  case  without  a  pseudo-magnetic  field  present;  i.e. 
the  existence  of  edge  modes  in  system  (4).  In  this  case  we  found  a  class  of  stationary  edge  cor¬ 
responding  to  so-called  zig-zag  boundary  conditions  cf.  [36]  and  references  therein.  In  [37]  we 
developed  an  asymptotic  theory  that  describes  and  extends  the  results  of  [11]  to  the  case  where 
there  is  a  pseudo-magnetic  field.  The  advantage  of  the  asymptotic  theory  is  that  all  formulae  are 
explicit  which  allows  us  to  readily  generalize  to  other  cases  of  physical  interest. 

In  addition  to  studying  the  dynamics  and  propagation  of  edge  wave  dynamics  in  honeycomb 
lattices  we  also  analyzed  the  propagation  of  multidimensional  nonlinear  surface  plasmon  waves 
in  dielectric/metal  interfaces  [38].  Knowledge  of  both  areas  is  important  in  order  for  us  to  deeply 
understand  the  nature  of  wave  dynamics  along  edges  and  surfaces. 

Using  the  HC  lattice  theory  with  a  pseudo-field  we  were  able  to  examine  situations  with  more 
general  pseudo-magnetic  fields.  We  find  that  the  linear  edge  modes  are  remarkably  robust  in  a 
certain  regime  of  parameter  space;  in  other  regimes  they  disperse.  We  have  also  used  the  theory  to 
generalize  to  nonlinear  lattices  and  have  found  that  the  envelope  function  associated  with  the  edge 
modes  satisfy  the  classical  nonlinear  Schrodinger  (NLS)  equation.  We  will  continue  studying  this 
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problem  and  will  carefully  investigate  whether  the  nonlinear  soliton  states  can  propagate  without 
dispersion  along  boundaries  even  when  there  are  corners  present. 

The  robust  behavior  of  the  propagating  edge  modes  appears  to  be  correlated  with  underlying 
topology  associated  with  the  Bruillouin  zone  [11,  39].  Propagating  edge  modes  persist  over  long 
distances  despite  the  fact  that  we  are  in  a  gap  regimes  of  the  edge  mode  spectrum.  This  situa¬ 
tion/effect  is  sometime  referred  to  as  a  ‘Floquet  topological  insulator’  in  the  physics  literature. 
Topological  insulators  are  insulators  in  the  bulk  but  admit  boundary  conducting  sates  whose  en¬ 
ergies  lie  in  a  bulk  energy  gap.  The  term  Floquet  derives  from  the  fact  that  the  pseudo-field  in 
periodic  in  the  longitudinal  direction.  Preliminary  calculations  indicate  that  the  robustness  of  the 
Floquet  topological  insulator  in  the  linear  case  can  also  lead  to  similar  results  for  the  nonlinear 
regime.  We  are  continuing  to  study  this  in  detail. 

PT-Symmetric  systems 

Another  area  of  interest  in  modern  optics  is  so  called  VT  systems.  These  systems  are  made 
by  introducing  gain  and  loss  in  just  such  a  way  that  the  equations  describing  the  propagation  of 
fields  which  are  invariant  under  the  combined  action  of  spatial  inversion,  V,  and  time  reversal, 
T.  'PT-symmctric  optical  systems  represent  a  new  class  of  optical  metamaterials  allowing  for 
greater  control  of  light.  In  particular,  VT  -symmetric  optical  systems,  phenomena  such  as  novel 
beam  diffraction  patterns  and  unidirectional  invisibility  have  been  observed.  All  of  this  phenomena 
depends  on  the  fundamental  question  of  whether  spectra  is  real  or  if  a  ‘phase  transition’  [40]  has 
caused  spectra  with  non-zero  imaginary  part  to  appear.  Establishing  whether  phase-transitions 
do  or  do  not  occur  has  been  studied  in  many  contexts  cf.  [41].  In  [42],  starting  from  a  linear 
Schrodinger  equation  of  the  type  discussed  above,  see  eq.  (1)  with  a  PT  symmetric  lattice  potential, 
we  develop  a  necessary  and  sufficient  condition  which  establishes  when  phase  transitions  can  or 
can  not  occur  in  VT  symmetric  honeycomb  lattices  for  small  VT  perturbations.  We  we  show  that 
honeycomb  potentials  with  added  symmetry  allow  us  to  find  VT  perturbations  which  satisfy  this 
condition.  Via  numerical  experiments,  it  is  seen  that  VT  -symmetric  lattices  satisfying  the  analytic 
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condition  do  not  exhibit  phase  transitions  for  a  range  of  parameter  regimes,  sometimes  even  for 
relatively  large  VT  perturbations.  This  goes  beyond  the  standard  theory,  and  shows  with  added 
symmetry  that  carefully  designed  VT  symmetric  lattices  can  be  robust  against  phase  transitions. 

In  a  related  direction  we  have  found  that  a  ‘simple  looking’  nonlocal  PT-symmtric  equation  is 
integrable.  The  equation  takes  the  form 

iqt(x,t )  =  qxx(x,t)  ±  2q(x,t)q*(~x,t)q(x,t)  (8) 

where  *  denotes  complex  conjugation  and  q(x,t )  is  a  complex  valued  function  of  the  real  vari¬ 
ables  x  and  t.  Eq.  (8)  admits  a  linear  (Lax)  pair  formulation  and  possesses  an  infinite  number  of 
conservation  laws;  hence  it  is  an  integrable  system.  Using  the  the  inverse  scattering  transform  we 
can  linearize  the  equation,  corresponding  to  rapidly  decaying  initial  data,  and  obtain  solutions  to 
Eq.  (8)  including  pure  solitons  solutions  [43].  Some  of  the  important  properties  of  the  nonlocal 
NLS  equation  are  contrasted  with  the  classical  NLS  equation  where  the  nonlocal  nonlinear  term 
q*(—x,  t)  is  replaced  by  q*  (x,  t ) .  Indeed  we  note  that  both  equation  (8)  and  the  classical  NLS  share 
the  symmetry  that  when  x  — »  —x,  t  — >•  —t  and  a  complex  conjugate  is  taken,  then  the  equation 
remains  invariant.  Thus,  the  new  nonlocal  equation  is  PT  symmetric  which,  in  the  case  of  classi¬ 
cal  optics,  amounts  to  the  invariance  of  the  induced  potential  V (x,  t)  =  q(x,  t)q*(—x,  t )  under  the 
combined  action  of  parity  and  time  reversal  symmetry.  Finally,  wave  propagation  in  PT  symmet¬ 
ric  coupled  waveguides/photonic  lattices  has  been  experimentally  observed  in  classical  optics  cf. 
[40], 

Similarly  we  find  the  following,  also  ‘simple  looking’,  nonlocal  discrete  nonlinear  Schrodinger 
equation  is  integrable  [44] 


/<l<(ljl  =  j ^  (?n+ 1  -  2 qn  +  qn- 1)  ±  qnq*-n  {qn+i  +  qn- 1) 


(9) 


where  *  denotes  complex  conjugation  and  qn  is  a  complex  function  and  n  is  an  integer.  As  with 
equation  (8)  this  discrete  equation  (9)  admits  a  linear  (Lax)  pair  formulation  and  possesses  an 
infinite  number  of  conservation  laws;  hence  it  too  is  an  integrable  system.  Corresponding  to  rapidly 
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decaying  initial  data  one  can  exactly  linearize  Eq.  (9)  using  the  inverse  scattering  transform  and 
obtain  solutions  including  pure  solitons.  Some  of  the  important  properties  of  Eq.  (9)  are  contrasted 
with  the  better  known  model  where  the  nonlocal  nonlinear  term  qn<f-n  is  replaced  by  \qn\2.  Indeed 
as  with  the  continuous  model,  the  discrete  equation  shares  the  symmetry  that  when  n  — *  —n,  t  — > 
—t  and  a  complex  conjugate  is  taken,  then  the  equation  remains  invariant.  Thus,  the  new  discrete 
nonlocal  equation  is  PT  symmetric  which,  in  the  case  of  classical  optics,  amounts  to  the  invariance 
of  the  self-induced  potential  Vn  =  qnq*-n  under  the  combined  action  of  parity  and  time  reversal 
symmetry. 
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Dynamics  of  ultra-short  mode-locked  laser  pulses 

Research  developments  with  mode-locked  lasers,  such  as  Ti: sapphire  lasers,  have  enabled  sci¬ 
entists  to  generate  regularly  spaced  trains  of  ultrashort  pulses,  which  are  separated  by  one  cavity 
round-trip  time.  Fig.  1 1  below  shows  a  schematic  of  a  mode-locked  Ti:sapphire  laser  and  the 
emitted  pulse  train. 


Figure  1 1:  Ti:sapphire  laser  (left)  and  the  emitted  pulse  train  (right). 


Associated  with  the  spectrum,  or  the  Fourier  transform  of  the  pulse  train,  is  a  frequency  comb 
whose  frequencies  are  separated  by  the  laser’s  repetition  frequency.  Progress  in  the  development 
of  optical  oscillators  has  been  made  possible  by  controlling  these  femtosecond  frequency  combs. 
Extremely  stable  frequency  combs  have  been  generated  by  Ti:s  laser  systems,  but  other  types  of 
mode-locked  lasers  such  as  SnForsterite  and  fiber  lasers  are  also  being  actively  studied.  It  is 
important  to  have  useful  mathematical  models  of  these  laser  systems. 

In  our  research  we  have  been  studying  a  distributed  dispersion-managed  equation  which  we 
term  the  power  energy  saturation  (PES)  model.  For  a  pulse  amplitude  u(z ,  £),  power  P(z,  t)  =  |u|2, 
and  energy  E(z)  =  |m|2  df,  propagating  in  the  2  direction,  the  normalized  or  dimensionless 

equation  we  study  takes  the  form 


.du  d(z)d2u 
l!h+  2  dt2 


+  n{z)\u\2u 


ig  ir  il 

1  +  E/Esa.tU  +  1  +  E/EsatUtt  ~  1  +  P/PsatU 


(10) 
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where  the  constant  parameters  g,  r,  /,  Esat,  Psat  are  positive.  The  first  term  on  the  right  hand 
side  represents  saturable  gain,  the  second  is  nonlinear  filtering  (r  ^  0)  and  the  third  is  saturable 
loss.  This  model  generalizes  the  well-known  master  laser  equation  originally  developed  by  Haus 
and  collaborators.  When  the  loss  term  is  approximated  in  the  weakly-nonlinear  regime  by  a  first 
order  Taylor  polynomial  we  obtain  the  master  laser  equation.  Hence,  the  master  laser  equation  is 
included  in  the  power  saturated  model  as  a  first  order  approximation. 

In  our  papers  we  discuss  in  detail  how  this  equation  applies  to  the  system  described  in  Fig.  11. 
The  system  is  dispassion  managed  because  there  is  a  difference  between  the  dispersion  inside  the 
crystal  and  in  the  prisms/mirrors. 

In  our  earlier  research  investigations  in  fiber  optics  [45]  we  derived,  based  on  the  asymptotic 
procedure  of  multiple  scales,  a  nonlinear  integro-differential  equation  (not  given  here  due  to  space 
considerations)  which  governs  the  dynamics  of  dispersion-managed  pulse  propagation.  This  gov¬ 
erning  equation  is  referred  to  as  the  dispersion-managed  nonlinear  Schrodinger  (DMNLS)  equa¬ 
tion.  When  there  is  no  gain  or  loss  in  the  system,  for  strongly  dispersion-managed  systems,  the 
DMNLS  equation  plays  the  role  of  the  “pure”  NLS  equation-which  is  the  relevant  averaged  equa¬ 
tion  when  there  is  either  small  or  no  dispersion-management.  When  gain  and  loss  are  included  as 
in  the  PES  equation  we  find  a  modification  of  the  “pure”  DMNLS  equation  [46]. 

With  or  without  dispersion-management  the  PES  equation  naturally  describes  the  locking  and 
evolution  of  pulses  in  mode-locked  lasers  that  are  operating  in  the  soliton  regime.  In  our  research 
[46,  47]  we  took  a  unit  gaussian  input,  typical  values  of  the  parameters,  and  we  varied  the  gain 
parameter  g.  When  g  <  g*  no  localized  solution  was  obtained;  i.e.  in  this  case  the  effect  of  loss 
is  stronger  than  a  critical  gain  value  and  the  evolution  of  a  Gaussian  profile  decays  to  the  trivial 
solution.  Conversely,  when  g  >  g*,  there  exists  a  single  localized  solution,  u  =  U0(t)  exp(igz) 
where  fi,  called  the  propagation  constant,  is  uniquely  determined  given  the  specific  values  of  the 
other  parameters. 
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A  typical  situation  is  described  by  the  evolution  of  the  pulse  peak  for  different  values  of  the 
gain  parameter  g  as  shown  in  Fig.  12;  here  we  keep  all  terms  constant  and  only  change  the  gain 
parameter  g.  When  g  =  0.1  the  pulse  vanishes  quickly  due  to  excessive  loss  with  no  noticeable  os- 


Figure  12:  Evolution  of  the  pulse  peak  of  an  arbitrary  initial  profile  under  PES  with  different  values 
of  gain.  The  damped  pulse-peak  evolution  is  shown  with  a  dashed  line.  In  the  right  hand  figure, 
the  complete  evolution  is  given  for  g  =  0.3. 


dilatory  behavior;  the  pulse  simply  decays,  resulting  in  damped  evolution.  When  g  =  0.2, 0.3,  due 
to  the  loss  in  the  system  the  pulse  initially  undergoes  a  relative  to  its  amplitude  a  modest  decrease. 
However,  it  recovers  and  evolves  into  stable  solution.  Interestingly,  e.g.  when  g  =  0.7, 1,  and 
the  perturbations  are  not  small,  a  stable  evolution  is  nevertheless  again  obtained,  although  some¬ 
what  different  from  the  case  above.  Now  with  excessive  gain  in  the  system,  the  pulse  amplitude 
increases  develops  oscillations  but  a  steady  state  is  rapidly  reached. 

Generally  speaking  in  the  PES  model,  the  mode-locking  effect  is  present  for  g  >  g*,  a  critical 
gain  value.  Without  enough  gain  i.e.  g  <  g*,  pulses  dissipate  to  the  trivial  zero  state.  Furthermore, 
we  do  not  find  complex  radiation  states  or  states  whose  amplitudes  grow  without  bound  for  pa¬ 
rameter  regimes  we  studied.  In  terms  of  solutions,  Eq.  (10)  admits  soliton  states  for  all  values  of 
g  >  g*  >  l  (recall,  here  l  =  0.1).  This  was  also  shown  to  be  the  case  when  we  employed  analytical 
methods  (soliton  perturbation  theory). 

It  is  useful  to  note  that  only  for  a  narrow  range  of  parameters  does  the  master  laser  equation 
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(when  the  loss  term  is  taken  to  be  the  first  two  terms  of  the  Taylor  expansion  of  the  last,  power 
saturated  term,  in  the  PES  equation)  have  stable  soliton  solutions  or  mode-locking  evolution.  In 
general  the  solitons  are  found  to  be  unstable;  either  dispersing  to  radiation  or  evolving  into  nonlo- 
calized  quasi-periodic  states.  For  different  parameters,  the  amplitude  can  also  grow  rapidly  under 
evolution.  Thus,  the  basic  master  laser  equation  captures  some  qualitative  aspects  of  pulse  propa¬ 
gation  in  a  laser  cavity;  however,  since  there  is  only  a  small  range  of  the  parameter  space  for  which 
stable  mode-locked  soliton  pulses  exist,  it  does  not  reflect  the  wide  ranges  of  operating  conditions 
where  mode-locking  occurs. 

Interestingly,  as  the  gain  becomes  stronger  additional  soliton  states  are  possible  and  2,  3,  4  or 
more  coupled  pulses  are  found  to  be  supported.  This  means  that  strings  of  soliton  states  can  be 
obtained  [48]. 

Power  saturation  models  also  arise  in  other  problems  in  nonlinear  optics  and  are  central  in  the 
underlying  theory.  For  example,  power  saturation  models  are  important  in  the  study  of  the  dy¬ 
namics  of  localized  lattice  modes  (solitons,  vortices,  etc)  propagating  in  photorefractive  nonlinear 
crystals.  If  the  nonlinear  term  in  these  equations  was  simply  a  cubic  nonlinearity,  without  satu¬ 
ration,  two  dimensional  fundamental  lattice  solitons  would  be  vulnerable  to  blow  up  singularity 
formation,  which  is  not  observed.  Thus  saturable  terms  are  crucial  in  these  problems. 

In  recent  work  we  analyzed  how  the  mode-locking  mechanism  responds  to  higher  order  per¬ 
turbative  terms  including  third  order  and  nonlinear  dispersion  and  Raman  gain;  we  considered  the 
following  perturbation  of  the  PES  equation  (10)  with  constant  dispersion  and  nonlinearity 
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where  f3  represents  third  order  dispersion,  7  nonlinear  dispersion  or  self- steepening  and  R  the 
Raman  gain.  When  (3  =  7  =  R  =  0  this  equation  is  the  PES  model.  In  [49]  mode-lockiing 
was  found  to  be  maintained  even  with  these  higher  order  terms  terms;  mode-locking  occurs  for 
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both  the  anomalous  and  normal  regimes.  In  the  anomalous  regime,  these  perturbations  are  found 
to  affect  the  speed  but  not  the  structure  of  the  locked  pulses.  The  pulses  behave  like  solitons  of 
a  classical  nonlinear  Schrodinger  equation  and  as  such  a  soliton  perturbation  theory  verifies  the 
numerical  observations.  In  the  normal  regime,  the  effect  of  the  perturbations  is  small,  in  line  with 
recent  experimental  observations  [50]. 
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Dark  Solitons  and  mode  locked  lasers 

Dark  solitons  or  envelope  solitons  having  the  form  of  density  dips  with  a  phase  jump  across 
their  density  minimum,  are  fundamental  nonlinear  excitations  which  arise  for  example,  in  the 
defocusing  nonlinear  Schrodinger  (NLS)  equation.  They  are  termed  black  if  the  density  minimum 
is  zero  or  grey  otherwise.  The  discovery  of  these  structures,  which  dates  back  to  early  1970s, 
was  followed  by  intensive  study  both  in  theory  and  in  experiment:  in  fact,  the  emergence  of  dark 
solitons  on  a  modulationally  stable  background  is  a  fundamental  phenomenon  arising  in  diverse 
physical  settings.  Dark  solitons  have  been  observed  and  studied  in  numerous  physical  contexts 
including:  discrete  mechanical  systems  electrical  lattices,  magnetic  films,  plasmas,  fluids,  atomic 
Bose-Einstein  condensates  as  well  as  nonlinear  optics. 

In  nonlinear  optics,  dark  solitons  have  some  advantages  as  compared  to  their  bright  coun¬ 
terparts,  which  are  solutions  of  the  focusing  NLS  mode.  Dark  solitons  are  more  robust  against 
Gordon-Haus  jitter,  background  noise  and  higher  order  dispersion. 

Recently,  there  has  been  an  interest  in  “dark  pulse  lasers”,  namely  laser  systems  emitting  trains 
of  dark  solitons  as  envelopes  of  a  continuous  wave  (cw)  emitted  by  the  laser;  various  experimental 
results  have  been  reported  utilizing  fiber  ring  lasers,  quantum  dot  diode  lasers  and  dual  Brillouin 
fiber  lasers  [51,  52].  These  works,  apart  from  introducing  a  method  for  a  systematic  and  control¬ 
lable  generation  of  dark  solitons,  can  potentially  lead  to  other  important  applications  related,  e.g. 
with  optical  frequency  combs,  optical  atomic  clocks,  fiber  optics  etc.  An  important  aspect  in  these 
studies  is  the  ability  of  the  laser  system  to  induce  a  fixed  phase  relationship  between  the  modes  of 
the  laser’s  resonant  cavity,  i.e.  to  mode-lock ;  in  such  a  case,  interference  between  the  laser  modes 
in  the  normal  dispersion  regime  causes  the  formation  of  a  sequence  of  dark  pulses  on  top  of  the 
stable  cw  background  emitted  by  the  laser. 

We  have  studied  dark  solitons  subject  to  general  perturbations 

iuz  -  ^ utt  +  \u\2u  =  F[u ]  (12) 

where  \F\  <C  1.  We  have  developed  a  general  perturbation  theory  for  dark  solitons  [53].  One  of 
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our  results  is  that,  in  general,  perturbations  of  dark  solitons  induce  a  small  but  wide  shelf.  We  can 
calculate  the  perturbation  to  the  soliton  parameters  as  well  as  the  shelf. 

As  a  prototypical  situation  involving  dark  solitons,  we  considered  the  perturbation  theory 
within  the  framework  of  the  PES  equation  (10),  in  the  normal  regime  [54],  We  found  that  general 
initial  conditions  evolve/mode-lock  into  dark  solitons  under  appropriate  gain  requirements.  The 
resulting  pulses  are  essentially  dark  solitons  of  the  unperturbed  nonlinear  Schrodinger  equation. 
In  a  ring  like  laser  system  the  PES  model  also  develops  a  rather  pronounced  shelf.  In  Figure  (13) 
below,  a  typical  situation  is  depicted  for  a  ring  laser  system  described  by  the  PES  equation.  The 
shelf  is  seen  to  be  prominent  in  the  PES  equation.  No  shelf  exists  in  the  unperturbed  NLS  equation. 
In  the  inset  of  the  figure  the  phase  change  across  the  dark  soliton  is  given  for  the  PES  model  vs.  the 
unperturbed  NLS  equation.  In  practice  the  shelf  might  interact  with  background  noise  creating  a 
complex  background  upon  which  the  dark  soliton  resides.  Nevertheless  dark  solitons  are  extremely 
stable  even  with  such  background  disturbances. 
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Figure  13:  The  development  of  a  shelf  in  the  solution  to  the  PES  equation  in  a  ring  laser  config¬ 
uration  (solid  line)  as  compared  to  the  solution  of  the  NLS  equation  (dashed  line).  The  shelf  is 
evident  in  the  amplitude.  The  inset  shows  how  the  phase  changes  across  the  soliton. 
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Dispersive  Shock  Waves 

Shock  waves  in  compressible  fluids  is  a  classically  important  field  in  applied  mathematics  and 
physics,  whose  origins  date  back  to  the  work  of  Riemann.  Such  shock  waves,  which  we  refer  to  as 
classical  or  viscous  shock  waves  (VSWs),  are  characterized  by  a  localized  steep  gradient  in  fluid 
properties  across  the  shock  front.  Without  viscosity  one  has  a  mathematical  discontinuity;  when 
viscosity  is  added  to  the  equations,  the  discontinuity  is  “regularized”  and  the  solution  is  smooth. 
An  equation  that  models  classical  shock  wave  phenomena  is  the  Burgers  equation 

ut  +  uux  =  uuxx  (13) 

If  v  =  0,  we  have  the  inviscid  Burgers  equation  which  admits  wave  breaking.  When  the  underlying 
characteristics  cross  a  discontinuous  solution,  i.e.  a  shock  wave,  is  introduced  which  satisfies 
the  Rankine-Hugoniot  jump  conditions  which,  in  turn,  determines  the  shock  speed.  Analysis  of 
Burgers  equation  shows  that  there  is  a  smooth  solution  given  by 


1  1  ,  f  1  ,  1  V, 

u  — - tanhl  —lx - 1)  1 

2  2  l4z/  2  n 

which  tends  to  the  shock  solution  as  v  — >■  0.  Thus  the  mathematical  discontinuity  is  regularized 
when  viscosity  v  is  introduced.  A  typical  regularized  shock  wave  can  be  seen  in  Fig.  14  below  in 
red. 


Another  type  of  shock  wave  is  a  so-called  dispersive  shock  wave  (DSW).  Early  observations  of 
DSWs  were  ion-acoustic  waves  in  plasma  physics;  indeed  the  Korteweg-deVries  (KdV)  equation 
describes  ion-acoustic  waves  in  plasmas.  Subsequently,  Gurevich  and  Pitaevskii  [55]  studied  the 
small  dispersion  limit  of  the  KdV  equation.  They  obtained  an  analytical  representation  of  a  DSW. 
As  opposed  to  a  localized  shock  as  in  the  viscous  problem,  the  description  of  a  DSW  is  one  with  a 
sharp  front  with  an  expanding,  rapidly  oscillating  rear  tail.  The  Korteweg-deVries  (KdV)  equation 
with  small  dispersion  is  given  by 


tif  +  uux  =  eu3 


(14) 
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where  |e|  -C  1  regularizes  the  discontinuity  that  otherwise  would  be  present.  The  mathematical 
technique  used  to  analyze  DSWs  relies  on  wave  averaging,  often  referred  to  as  Whitham  theory. 
Whitham  theory  is  used  to  construct  equations  for  the  parameters  associated  with  slowly  varying 
wavetrains;  it  provides  an  analytical  basis  for  DSW  dynamics.  For  KdV  the  Whitham  equations 
can  be  transformed  into  Riemann-invariant  form,  which  can  be  analyzed  in  detail.  The  picture  and 
details  are  quite  different  from  viscous  shock  waves. 

A  viscous  shock  wave  which  occurs  in  Burgers  equation  is  depicted  below  in  Fig.  14-in  red; 
a  typical  DSW  associated  with  the  KdV  equation  is  illustrated  in  blue.  For  the  KdV  DSW  one 
finds  that  there  are  two  speeds  associated  with  a  DSW:  one  is  the  speed  associated  with  the  frontal 
wave  which  is  a  soliton  (located  at  xs  in  the  figure),  and  the  other  speed  corresponds  to  the  group 
velocity  of  near  linear  trailing  waves  on  the  rear  end  (depicted  by  xt  in  the  figure)  of  the  DSW.  This 
is  very  different  from  the  classical  or  viscous  shock  wave  (located  at  xc  in  the  figure)  associated 
with  the  Burgers  equation.  Interestingly,  the  structure  of  the  KdV  DSW  is  strikingly  similar  to  the 
original  plasma  observations. 


Figure  14:  Left  figure:  typical  DSW  satisfying  the  KdV  eq.  (14)  and  a  classical  shock  wave 
satisfying  Burgers  eq.  (13);  right  figure:  typical  blast  wave-experiment;  numerical  simulation  is 
given  below  the  experiment. 


Recent  experiments  in  Bose-Einstein  condensates  (BEC)  and  nonlinear  optics  have  enhanced 
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interest  in  DSWs.  The  BEC  experiments,  originally  performed  in  the  Physics  Department  at  the 
University  of  Colorado,  motivated  our  studies  [56].  Recent  experiments  in  nonlinear  optics  carried 
out  in  the  laboratory  of  J.  Fleischer  at  Princeton  University  have  also  observed  similar  blast  waves 
and  other  interesting  DSW  phenomena  [57].  Since  the  equations  governing  the  phenomena  are 
similar,  this  provides  further  support  for  our  DSW  theory. 

While  interactions  of  viscous  shock  waves  are  well  known,  the  understanding  of  interacting 
DSWs  is  still  at  an  early  stage.  Nevertheless  we  have  made  some  progress  towards  understanding 
the  interaction  picture.  We  note  that  in  recent  nonlinear  optics  experiments  [57]  interacting  DSWs 
were  observed.  We  are  investigating  DSW  interactions  in  physically  interesting  systems  by  em¬ 
ploying  both  Whitham  methods  and  asymptotic  analysis  applied  to  the  solution  obtained  by  the 
inverse  scattering  transform. 

Recently  we  have  investigated  the  KdV  equation  with  step-like  data.  We  found  that  the  long- 
time-asymptotic  solution  of  the  KdV  equation  for  general,  step-like  data  is  a  single-phase  DSW; 
this  DSW  is  the  largest  possible  DSW  based  on  the  boundary  data.  We  find  this  asymptotic  solu¬ 
tion  using  the  inverse  scattering  transform  and  matched  asymptotic  expansions  and  confirmed  the 
phenomena  numerically.  So  while  multi-step  data  evolve  to  have  multiphase  dynamics  at  interme¬ 
diate  times,  these  interacting  DSWs  eventually  merge  to  form  a  single-phase  DSW  at  large  time 
[58,59,  60], 

Dispersive  shock  waves  are  an  interesting  and  developing  area  of  research  which  we  believe 
will  play  an  increasingly  important  role  in  nonlinear  optics  applications  and  other  areas  of  physics. 
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7.  South  African  Numerical  and  Applied  Mathematics  Symposium,  Stellenbosch, 
SA,  Nonlinear  waves  from  beaches  to  dispersive  shock  waves,  April  3,  2013 

8.  Conference  Board  on  Math  Sci.  conference  on  two-dimensional  solitons,  water 
waves;  Edinberg,  Texas,  May  20-24,  2013;  Two-dimensional  water  waves  theory 
and  ocean  observations,  May  20,  2013 

9.  BIRS  Workshop  on  Water  Waves:  Banff  Canada,  June  30-July  5,  2013;  Two- 
dimensional  water  waves  theory  and  ocean  observations.,  July  1,  2013 

10.  The  8th  Symposium  on  Quantum  Theory  and  Symmetries,  Mexico  City  Aug.  5-9, 

2013,  Nonlinear  waves  from  beaches  to  dispersive  shock  waves,  Aug.  5,  2013 

1 1 .  Deformations,  Localized  Edge  States  and  Associated  Pulse  Dynamics  in  Nonlinear 
Optical  Honeycomb  Lattices,  AFOSR  Workshop:  Nonlinear  Optics,  Arlington, 
VA 

12.  Electrical  Engineering  Dept.  Colloquium,  North  Carolina  State  Univ,  Jan.  17, 

2014,  “Nonlinear  waves  from  beaches  to  photonic  lattices” 

13.  Applied  Mathematics  Dept.,  Colloquium,  Calif.  Inst,  of  Tech.,  Feb  3,  2014,  ‘Non¬ 
linear  waves  from  beaches  to  photonic  lattices” 
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14.  Center  for  Nonlinear  Studies,  Los  Alamos  National  Laboratory,  Colloquium,  March 
10,  2014,  ‘Nonlinear  waves  from  beaches  to  photonic  lattices” 

15.  Frontiers  in  Applied  Mathematics,  New  Jersey  Institute  of  Technology,  May  22- 
23,  2014,  May  22,  2014  “Nonlinear  waves,  solitons  and  applications” 

16.  Nonlinear  Waves  in  Sicily,  June  8-13,  2014,  June  9,  2014  “Nonlinear  wave  in 
photonic  lattices” 

17.  SIAM  Conference  NL  Waves  and  Coherent  Structures,  Cambridge  University,  UK, 
Aug.  11-14,  2014,  Aug.  11,  2014,  “Nonlinear  waves  from  beaches  to  photonics” 

18.  Nonlinear  Optics  Workshop,  Air  Force  BRIC,  Arlington,  VA.  Oct.  2,  2014,  “Pho¬ 
tonic  honeycomb  lattices:  localized  linear  and  nonlinear  edge  states” 

19.  University  of  Ioannina,  “Mathematics  and  International  Collaborations”,  Oct.  22, 
2014 

20.  University  of  Ioannina,  “Nonlinear  Waves  -  Still  Many  Surprises”,  Oct.  23,  214 

21.  University  of  Colorado,  Colorado  Springs,  Colorado  Nonlinear  days,  “Nonlinear 
waves  -  Always  Alluring”,  Nov.  1,  2014 

-  Consultative  and  Advisory  Functions  to  Other  Laboratories  and  Agencies:  None 
NEW  DISCOVERIES,  INVENTIONS,  OR  PATENT  DISCLOSURES:  None 

HONORS/AWARDS: 

-  Named  as  one  of  the  most  highly  cited  people  in  the  field  of  Mathematics  by  the  ISI 

Web  of  Science,  2003-present 

-  Named  Fellow  of  the  American  Mathematical  Society  (AMS)  2012- 

-  Martin  Kruskal  Prize/Lecture:  Aug.  11,  2014,  Cambridge  University,  Awarded  by 

SIAM,  Activity  Group  on  Nonlinear  Waves  and  Coherent  Structures 

-  Doctor  Honoris  Causa,  Oct.  2014,  University  of  Ioannina,  Greece 
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